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We study here a standard next-nearest-neighbor (NNN) model of ballistic growth on one- and 
two-dimensional substrates focusing our analysis on the probability distribution function P(M, L) 
. of the number M of maximal points (i.e., local "peaks") of growing surfaces. Our analysis is based 

on two central results: (i) the proof (presented here) of the fact that uniform one-dimensional 
ballistic growth process in the steady state can be mapped onto "rise- and- descent" sequences in the 
ensemble of random permutation matrices; and (ii) the fact, established in Ref. [22ll . that different 
«4 ■ characteristics of "rise-and-descent" patterns in random permutations can be interpreted in terms 

of a certain continuous-space Hammersley-type process. For one-dimensional system we compute 
P(M, L) exactly and also present explicit results for the correlation function characterizing the 
enveloping surface. For surfaces grown on 2d substrates, we pursue similar approach considering 
the ensemble of permutation matrices with long-ranged correlations. Determining exactly the first 
three cumulants of the corresponding distribution function, we define it in the scaling limit using an 
. expansion in the Edgeworth series, and show that it converges to a Gaussian function as L — > oo. 

e 

-6 

C ■ 

Within the recent years much effort has been devoted to theoretical analysis of properties of surfaces obtained 
• by aggregation of particles. Several models describing various properties of clusters grown by different deposition 
processes have been proposed. To name but a few, we mention the famous Kardar-Parisi-Zhang (KPZ) pj and the 
Edwards- Wilkinson (EW) models Q , models of surfaces grown by Molecular Beam Epi taxy (M BE) (see, for example, 
Ref. 4]), Polynuclear Growth (PNG) H H Q, @ and by Ballistic Deposition (BD) [jjj HI El El , in which case 
, particles are sequentially added to a growing surface along ballistic trajectories with random initial positions and 
qq ' specified direction. 



I. INTRODUCTION 



For these models of surface growth, a number of important theoretical achievements concerning the statistics of 



extrema has been made. In particular, in Ref. :3j the distributions of the maximal heights of the ID Edwards-Wilkinson 
' and of the KPZ interfaces were determined exactly. In Ref.Q it was realized that the height distribution of the PNG 
surfaces coincides with the so-called Tracy- Widom distribution [|J, which appears in the theory of random matrices. 
In Ref. EH it has been shown that in the thermodynamic limit BD exhibits the KPZ scaling behavior. Moreover, 
a discrete BD model has been shown recently to be a very convenient tool for studying non-Abelian entanglement 
properties of braided directed random walks 16]. Finally, it has been found that in many models of ballistic growth, 
as well as in their continuum-space counterparts belonging to the KPZ universality class the average velocity of 
cluster's growth is governed by the density of local minima of the enveloping surface |17|. 



C . 

In this work we analyze the structure of the enveloping surface in a standard next-nearest-neighbor (NNN) model of 
ballistic growth on one- and two-dimensional substrates. Interpreting this model in terms of permutations of the set 
1, 2, 3, . . . , L, where L is the number of lattice sites and the numbers drawn at random from this set determine local 
heights of the surface, we calculate the Probability Distribution Function (PDF) of the number of maximal points 
(i.e., local "peaks"). Our analysis is based on two central results: 



(i) A proof, presented in this paper, of the fact that a ballistic growth process in the steady state can be formulated 
exactly in terms of a "rise-and-descent" pattern in the ensemble of random permutation matrices; 
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(ii) An observation made in Ref. 22] that the "rise-and-descent" patterns can be treated efficiently using a recently 
proposed algorithm of a Permutation Generated Random Walk. 

We remark that the expected value and the variance of the number of local peaks in surfaces grown by ballistic 
deposition on a one-dimensional and on a two-dimensional honeycomb lattices have been already calculated by J. 
Desbois in Ref.|l8l| using a certain decoupling of the hierarchy of coupled differential equations describing evolution 
of the moments of higher order. This method provides correct results for the first two moments of the distribution 
function and, apparently, may be extended further for the calculation of higher moments by truncation of the higher 
level correlations. However, this approach is not completely rigorous. Our approach, on contrary, is mathematically 
exact, does not rely on any uncontrollable assumption and enables us to go beyond the results obtained in [l^. In 
particular, in one dimension we calculate a complete distribution function of the number of local peaks exactly. Apart 
of that, we also present explicit results for the correlation functions characterizing the enveloping surface. For surfaces 
emerging in two-dimensional ballistic growth, we reduce the problem to the analysis of the ensemble of permutation 
matrices with long-ranged correlations. Determining exactly the first three cumulants of the corresponding PDF, we 
obtain the distribution function in the scaling limit using expansions in the Edgeworth series. 

The paper is outlined as follows. In Section II we formulate our models in one- and two-dimensions, discuss, 
on an intuitive level, a relation between a sequential growth of patterns in ballistic aggregation and " dynamics" on 
permutations, and finally present the main results of this work. Section III is devoted to a rigorous description of 
the relation between ballistic deposition and an "updating" dynamics on permutations. Further on, in Section IV we 
analyze the probability distribution function of the number of local peaks on one-dimensional substrates. Here, we 
briefly outline the well-known results of Stembridge on the peak numbers [l^, describe the model and basic results 
obtained for the so-called Permutation Generated Random Walks and show how the operator formalism developed 
in this work can be extended for the calculation of the moments of the probability distribution function. Next, in 
Section V we present a derivation of an explicit expression of the probability of having two local surface peaks at 
distance I apart of each other (such that there is no peaks inbetween). In Section VI we focus on two-dimensional 
situation and show how our previous analysis can be extended to this case. We calculate exactly first three cumulants 
of the corresponding probability distribution function of having M peaks on a square lattice containing L sites and 
then, using an expansion in the Edgeworth series, show that this function converges to a Gaussian as L — > oo. Finally, 
in Section VII we conclude with a brief summary of our results. 



A standard one-dimensional ballistic deposition model with next-nearest-neighboring (NNN) interactions is formu- 
lated as follows (for more details, see Ref. 14]). Consider a box divided in L columns (of unit width each) enumerated 
by index i (i = 1,2, ....L). For simplicity, we assume the periodic boundary conditions, such that the leftmost and 
the rightmost columns are neighbors. 

At the initial time moment, (n = 0), the system is deemed empty. Then, at each tick of the clock, n = 1,2, N, 
we deposit an elementary cell ("particle") of unit height and width in a randomly chosen column. Suppose that the 
distribution on the set of columns is uniform. Define the height, h(i, n), in the column i at time moment n. Assume 
now, as it is depicted in Fig^ that the cells in the nearest neighboring columns interact in such a way that they can 
only touch each other by corners, but never by their vertical sides. This implies that after having deposited a particle 
to the column i, the height of this column is modified according to the following rule: 



If at the time moment n nothing is added to the column i, its height remains unchanged: h(i,n + 1) = h(i,n). A 
set of deposited particles forms a pile as shown in Fig^i for L — 6 columns and N — 6 particles. Here, for example, 
h(l,6) = 1, h(2,6) =2, etc. 

Now, we call the local maxima of the pile as the "peaks" . More specifically, take the set H of heights at some time 
moment n: TL — {h(l, n), /i(2, n), h(L,n)}. We say that the column i contains a "peak" at time moment n if the 
height of this column satisfies the following two-sided inequality: 



II. MODELS, DEFINITIONS AND MAIN RESULTS 



A. Surfaces in standard ballistic growth process 



h(i, n + 1) = max{h(i — 1, n), h(i, n), h(i + 1, n)} + 1. 



(1) 



h(i, n) > h(i — 1, n) 
h(i, n) > h(i + 1, n). 



(2) 
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FIG. 1: Sequential growth of a heap in the ballistic deposition process. 



Note that in Fig^ there are two peaks situated in the columns i = 3 and i = 5. The collection of peaks T is the 
subset of Ti and forms the "roof" — the set of upmost (or "removable" 28]) particles. In Fig|T] peaks are denoted by 
gray squares and other particles - by white ones. The relation of this ballistic deposition process with the "updating 
dynamics" on permutations is schematically depicted in FigEl We briefly describe it below in Subsection B and in 
the Section III in more detail. 
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FIG. 2: An "updating dynamics" on permutations. 



In a similar fashion, the process of a two-dimensional ballistic deposition in a box with a square base L = m x m, 
(m is an integer), can be viewed as a sequential adding of elementary cubes in the columns satisfying the following 
rules (compare to |QJ): 

h(i,j,n + 1) = m&x{h(i - l,j,n), h(i + 1, j,n), h(i,j,n), h(i,j - l,n), h(i,j + l,n)} + 1, (3) 

where h(i,j,n) is the height of the column with coordinates (1 < < m) at deposition moment n (1 < n < 

2V). The cubes are added to the columns with the uniform distribution - see Fig|3] 

The "peak" of a two-dimensional landscape h(i, j, n) is defined as a local maximum in the set {h(i,j, n)} for some 
fixed moment n, if the following conditions are simultaneously fulfilled: 

' h(i,j,n) > h(i - l,j,n) 
h(i,j,n) > h(i + l,j,n) 

( 4 ) 

h(i,j,n) > h(i,j - l,n) 
, h(i,j,n) > h(i,j + l,n) 
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FIG. 3: Two-dimensional ballistic deposition and corresponding random landscape. 



Let us note that © - imply the periodic boundary conditions both in i and j coordinates (1 < {i, j} < to). The 
influence of the boundary condition on the expectation, variance and higher moments of peak numbers can be easily 
estimated and becomes negligible in the limit m — > oo. 

For these models, our goal is to evaluate the probability distribution function P(M, L) of having M peaks on a 
lattice containing L sites. As we proceed to show, in one-dimension this can be done exactly. Moreover, we are also 
able to calculate the "correlation function" p(l) defining the conditional probability that two peaks are separated by 
the interval I under the condition that the interval / does not contain peaks. In 2d, determining exactly first three 
cumulants of the PDF we define it in the scaling limit using expansion in the Edgeworth series and show that it 
converges to a Gaussian function as L — > oo. 



B. Interpretation of ballistic growth: "updating dynamics" and random permutations of natural series 

Let us emphasize that we are interested only in the statistics of peaks of a growing surface in the stationary state 
and disregard any questions concerning the statistics of heights. 

Dynamics of the set of peaks T in the ballistic deposition process depicted in Fig^ can be mapped onto the 
dynamics of "peaks" in the permutation matrix. We start by describing this connection on an intuitive level. To do 
this, let us proceed recursively. Suppose that we deposit a first particle in the column i\ of an L-column box. Next, 

ii-i 

take the row of L elements with " 1" at position i\ and "0" in all other places: (0, 0, 1, 0, 0). After dropping the 

» 2 -i 

second particle, say, in the column «2, take a row (0, ...,0, 1,0, ...,0) and place it over the first one creating a stack. 
Suppose that at some time n a new particle is added to the column which was occupied earlier, say, at time to, i.e. 
in = im {n > m)- It means that we have two identical rows in the stack. In this case, we remove the first of identical 
rows (i.e. deposited at time to) from the stack and eliminate the empty line by pulling down all rows deposited after 
time to, as it is depicted in Fig|2Jl. After some time, the stack will comprise L rows and, according to the described 
procedure, will not grow anymore but will be changed only due to updating of rows (by adding the new ones and 
by eliminating the old ones). By construction, this stack is an L x L permutation matrix. Connecting the nonzero 
elements in nearest-neighboring rows by a broken line, we can straightforwardly define the "descents" , "rises" and 
"peaks" in the permutation matrix — see FigEt- The number of peaks at time N in the L x L permutation matrix 
coincides then with the number of peaks in the heap after having deposited N particles in a box of L columns. The 
rigorous proof of this mapping is given in the Section IlIII 

Consequently, the described relation between the original ballistic deposition model and an "updating" dynamics 
on permutations, depicted in Fig|2i,d, allows us to view the characteristics of the surface obtained within the BD 
process from a different perspective: Consider, for example, a one-dimensional lattice containing L sites, on which we 
distribute numbers drawn randomly from the set 1, 2, 3, . . . , L. A number appearing at the site j determines the local 
"height" 29] of the surface. Now, we call this site as "a peak", if the number appearing on this site is bigger than the 
numbers on two neighboring sites — see Fig^l Generalization to two-dimensional square lattice with L = m x to, (to 
is an integer) sites is straightforward — see FigQJ): the only difference here is that we call as local "surface peaks" 
such sites j the numbers at which are bigger than numbers appearing at four adjacent sites. We prove rigorously 
in what follows that the probability of having M peaks on a lattice containing L sites on which we place numbers 
randomly drawn from the set 1, 2, 3, . . . , L, and the distribution function P(M, L) characterizing the number of local 
surface extrema in the surfaces of aggregates grown within the BD process, are identic. We note finally that the 
peak statistics in the latter model with natural numbers 1,2,3, ... ,L and in a related model with a lattice on sites 
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of which one places randomly numbers uniformly distributed in [0,1] has very interesting and unexpected features, as 
was recently communicated to us by B.Derrida |l5j . 
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(a) (b) 

FIG. 4: (a) One-dimensional (L = 6) and (b) two-dimensional square (1 = 4x4) lattices with periodic boundary conditions 
on sites of which we distribute numbers drawn from the set 1,2,3, ... ,L. These numbers determine the local heights of the 
surface. The sites with numbers bigger than those appearing on the neighboring sites are refereed to as the "peaks" . 



C. Main Results 



Main results of this work are as follows: 

(i) We prove that a uniform one-dimensional ballistic growth process in the steady state can be mapped onto the 
" rise-and-descent" sequences in the ensemble of random permutation matrices. 

(ii) We show that in one dimension the probability distribution function P(M, L) can be calculated exactly through 
its generating function 

oo 

W(s,L) = L\Y,s M+1 P(M,L), (5) 

M=0 

where 



M=0 



with A{M,L) being the so-called Eulerian numbers counting the number of "rises" (or "descents") in the 
ensemble of equally weighted permutation matrices of size L: 



<L + 1 



r=0 



A(M,£) = £(-lW ^ j(M-r) L , (7) 



where (?) denotes the binomial coefficient; here and henceforth we adopt a standard notation when (?) = b u^L 



(a-b)\ 



for b < a and (?) = otherwise. 



Inverting JSJ, we find the following exact expression for P(M, L): 



oL+2 M / L+l \ 1 m L + 1 

P(M, L) - E ( , _ + m) , (8) 

In the limit L —> oo, the PDF P(M, L) converges to a Gaussian distribution: 

3 IT f 45(M- ii) 2 ! 

Note that a similar result has been previously found in Ref.[22( for the sum of peaks and throughs in a per- 
mutation of length L within the context of the number of the U-turns made by the "Permutation Generated 
Random Walk" (PGRW). 



6 



(iii) We show that in one-dimensional systems the probability p(l) of having two peaks separated by the interval I 
under the condition that this interval does not contain other peaks, obeys: 



P(0 = 2 



,(* -!)(* + 2) 
(I + 3)! 



Curiously enough, this expression coincides with the distribution function of the distance between two "weak" 
bonds obtained by Derrida and Gardner [2^] in their analysis of the number of metastable states in a one- 
dimensional random Ising chain at zero's temperature. 

(iv) Using the cumulant expansion, we show that in two dimensions, in the limit L — > oo, the normalized probability 
distribution function p(x, L), where x = (M — fif D )/(j, converges to a Gaussian distribution 
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and first three cumulants are given by 
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(v) We have realized that local peaks in one- and two-dimensional equilibrium NNN ballistic deposition model 
are correlated and the correlation length extends over two lattice spacings: while by definition, the probability 
of having two peaks at the neighboring sites is zero, appearance of two peaks on next-nearest-neighboring 
sites (also on a diagonal for 2d) is higher than squared mean density of peaks (1/9 and 1/25 in Id and 2d, 
respectively). On contrary, the probability of having two peaks at a distance exceeding two lattice spacings is 
equal to the squared mean density, which signifies that they are uncorrelated. In other words, there are effective 
short-range "attractive" interactions between peaks. The corresponding weights of a few typical configurations 
are summarized in figure Fig0 



III. RIGOROUS DESCRIPTION OF THE RELATION BETWEEN "UPDATING DYNAMICS" ON 

PERMUTATIONS AND BD. 



As we have already mentioned in the previous section, dynamics of the set of peaks T in the ballistic deposition 
process depicted in Fig|21can be mapped onto the dynamics of "peaks" in the permutation matrix. Let us formulate 
this mapping explicitly. 

In a more rigorous approach, we have to consider two dynamical systems: one - on peak sets and the other - on 
permutations. Let us call the operation which corresponds to dropping a box in the column i on a peak set. That 
is if S is a peak set (i.e. a subset of {1, ... ,L} without two consecutive numbers) then 

Ei(S):=SU{i}/{i-l,i+l}. (10) 

The corresponding operation on permutations is defined by 

(cr(k) if a(k) < a{i) 
fi:=Fi(a) with fi(k):=iL if k = i (11) 

[cr(fc)-l if a(k) > a(i) 



One-Dimensional System 



Two-Dimensional System 
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P( • = 1/3 



P( 



)= 1/5 



P( ■ ■ M . . ) =0 



P( , 9 , 9 , ) = 2/15 > 1/3 x 1/3 



P( 



) = 



P( ■ * * ) = 1/3 x 1/3 



P( 



) = 1/25 



P( 



) = 2/45 > 1/25 



P( 



) = 1/20 > 1/25 



FIG. 5: Probabilities of having a single peak and a pair of isolated peaks some distance apart from each other in one- and 
two-dimensional models. Filled circles denote peaks - sites with numbers exceeding numbers appearing on nearest-neighboring 
sites. First column presents probabilities of several configurations in Id, while the second one describes the probabilities of a 
few possible configurations appearing in two dimensions. 



To proceed, we also need to define the reverse operation on permutations, Rk, which removes the upmost row and 
inserts it under the row k: 

(a(l) ifa(l)<k 

/i := R k (a) with fi(l) := Ik if a (I) = L (12) 

[o-(i) + 1 HL>a(l)>k 

It is clear now that Fi(cr) = \i if and only if Rk{n) = & where k = <r(i) is the row of the 1 in the column i. 
Denote now by Peak(cr) the peak set of the permutation a. Then, it is obvious that for all permutation a 

^(Peakfo-)) = Peak(Fi(cr)) for all i < L. (13) 

This simple but crucial observation allows us to translate the dynamics on permutations to the dynamics of peak 
sets. We will proceed in two steps: first, we will show that the limit probability measure of the dynamical system on 
permutations is equi-distributed and second, we will demonstrate that the image of this probability measure by the 
peak-set map is the limit probability measure on peak sets. 

Let (P(a)) a be a probability measure on permutations. Then, after one uniformly random chosen step Fi the new 
probability P'(fi) of a permutation fi is 

p '^) = iEE p ( ff ) ( 14 ) 

i— 1 a 

where the inner sum extends over the sets of all permutations a, such that Fi(a) — fi. However, by definition of Fi, 
if a permutation r is of the form t — Fi(a) then r(i) — L. Thus in the previous sum the only value of i such that 
there exists some a such that Fi(a) = fi is i = /i _1 (L), that is if i is the column where the 1 is in the top row L in 
fi. Now, here are exactly L permutations a verifying Fi{a) = /i, namely Rk{n) for k = 1 . . . L. Consequently, P'(/i) 
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is simply the average of the probabilities of the Rk(fi): 

P'(M) = iX>(iW)- ( 15 ) 
fc=i 

Moreover, it is clear that for all pairs of permutations a, \x there is a sequence (i r ) r of transformations which maps a 
on y«, that is fi = . . . F i3 F i2 F il (a). Consequently, the map F : P > P' is an irreducible Perron-Frobenius map whose 
maximal module eigenvalue is 1 with multiplicity 1. Thus, when the number of boxes tends to infinity, the probability 
converges to some limit and this limit is the unique, normalized (the sum of the coordinate is 1) eigenvector of F 
corresponding to the eigenvalue 1. In virtue of l|15(l . the uniform distribution on permutations is preserved by F so 
that it must be the limit distribution. 

Using it is possible to translate our permutations language to the language of the peak sets. Recall that we 
call a peak set any subset of {1, . . . , L} without two consecutive numbers. The only needed remark is that any non 
empty peak set S is the peak set of a permutation. Thus the map Peak extends to a surjective map from the set of 
probability measure on permutations to the set of probability measures on peak sets: 

Peak: Pi — ► Peak(P)(5) := ^ P{a) . (16) 

a : Pcak((r)=S 

Denote next E the map P(S) i— ► P'(S) where P'(S) is the new probability measure on sets after a uniformly random 
step Ei. Then lfH?|> can be written down as 

P(Peak(P)) = Peak(P(P)) , (17) 

for all probability measures on permutations P. In particular, the spectrum of E is included in the spectrum of F 
and moreover, the generalized eingenspaces (kernel of (F — \I) m for m large) of E are the image under Peak of the 
generalized eigenspaces of F. Consequently, E is a Perron-Frobenius map with maximal module eigenvalue 1 with 
multiplicity 1. As a consequence the limit distribution on peak sets exists and is the image by the map Peak of the 
limit distribution on permutations 

fllmit (S) = Peak (P limit (a)) (S) = i #{<r | Peak(a) = S} . (18) 
Before we proceed further, a few comments might be in order: 

First, as a matter of fact, this technique is quite general and can be applied to a more general notion of heaps of 
pieces. We demonstrate it on the following example. Let G = (V, E) be a finite oriented simple graph with vertex and 
edge set V and E and let L be the number of vertices. Each vertex v € V is associated with one "type" of pieces, and 
an arrow e = v i— > v' encodes the fact that when a piece of type v falls after a piece of type v', then it is placed over 
it, thus they are no more pieces of type v' at the top of the pile. The notion of permutation generalizes to the notion 
of standard labeling of the vertex of the graph, that is assignation of distinct numbers from 1 ... L to the vertices v. 
Then a vertex v is called a peak of a labeling I if l(v) > l(v') for all edges v t—> v' . One immediately notices that 
the previous reasoning applies, which implies that for all sets S of vertices the limit probability measure Pnmit('S') of 
S to be exactly the set of maximal pieces of a random heaps on the graph G is given by (|18(l . In this regard, the 
limit distribution of the 2D models can also be computed using these generalized peaks. Moreover, in the case of the 
oriented lines the notion of peaks reduces to the notion of a descents. 

Second, combinatorics of peaks of permutations has been recently reviewed by J. Stembridge in Ref.pjj- The peak 
algebra of Stembridge is a certain sub-algebra of the algebra of polynomials with infinitely many variables. It has a 
natural basis K$ indexed by peak sets S and therefore can serve as the support for generating series of probability of 
a peak set. In this regard, the generating series of the probability of each peak set has a very simple expression 

oo oo L 

J2t L ]T P(S)K s = J2jj E I Peak(er) = S}K S = exp(ifii) (19) 

L=0 SC{1...L} L=0 ' Sc{l...i} 

where K\ is the unique peak set where L = 1. Note that there is a minor difference between our work and Ref.[l9| 
since in the latter the extremities 1 and L were never considered as a peak, while we use the periodic boundary 
conditions. This difference should be unimportant for sufficiently large L. 

Furthermore it should be noticed that in pi) a different random walk is considered on peaks and conjectured on 
permutations with the same limit probability measure. Finally, the result of 20] suggests that this measure should be 
seen as some kind of generalized Plancherel measure associated with the degenerated Hecke-Clifford algebra instead 
of the symmetric group. 
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IV. PROBABILITY DISTRIBUTION FUNCTION P(M,L) IN ID 



We are now in position to determine exactly the probability P{M, L) that a random heap has a fixed number of 
peaks. From our previous analysis, it follows that such a probability is equal to the number B(M, L) of permutations 
of length L having exactly M peaks, divided by the total number of permutations L\, i.e. P(M,L) = B(M, V)j LI. 
Now, we recollect that if one considers permutation descents instead of peaks, the numbers A(M, L) of permutations 
of 1 ... L with exactly M descents are known as the Eulerian numbers, obey the following three-site recursion 



A(M, L) = (L- M + 1) A(M - 1, L - 1) + M A(M, L - 1) 
which leads to the following recurrence relation for the generating function U(t, L) = X)m=o t M+1 A(M, L) 



U(t,l)=t, U(t,L) 



t{l - t)jj-U(t, L - 1) + LtU(t, L-l) 



(20) 



(21) 



The recursion relation for peak-Eulerian number B(M, L) has been considered by Stembridge in [T]j, Remark 4.8. 
In our notations, his results attains the following form 

B(M, L) = (L - 2M + 2) B(M — 1,L— 1) + 2M B(M, L-l) (22) 

which in turn can be encoded by the recurrence relation obeyed by the generating series W{s, L) — 
Em=o sM+1 ^(^i) (note that compared to Stembridge there is no t k+1 term): 



4s 



W(s,l) = s, W{s,L) = 

Here is a table of a few first polynomials W(s,L): 

W(s,l 
W(s,2 
W(a,3 
W(s,4 
W(s,5 
W(s,6 
W{s,7 
W(s,8 
W(s,9 



2s(l 



LsW(s,L-l) 



(23) 



s 

2 s 
2s 2 - 
16 s 2 
16 s 3 



8 s z + 16 s 
272 s 3 + 416 s 2 + 32 s 
272 s 4 + 2880 s 3 + 1824 s 2 + 64 s 
7936 s 4 + 24576 s 3 + 7680 s 2 + 128 s 
7936 s 5 + 137216 s 4 + 185856 s 3 + 31616 s 2 



(24) 



256 s 



Differentiating polynomials W(s,L) and setting s = 1, we can find, in principle, any moments of the probability 
distribution function. For example, the expectation of the number of peaks is for L > 2; the variance is 2 ^ 2 for 
L > 4, while the third central moment is equal to — 2 g 4 1; 2 for L > 8. 



A. Permutation Generated Random Walks 



In this subsection we briefly outline the basic notions concerning the so-called Permutation Generated Random Walk 
(PGRW) proposed in Ref . |22j] . Consider a given permutation it = {tti, 7r 2 , 7r 3 , ttl+i} of [L + 1] = 1, 2, 3, . . . , L + 1 
and rewrite it as a 2-line table: 

/ 1 2 3 ... L+iy 

\ Tl 1"2 7T3 ■•• ^L+l / 

Suppose that this table assigns some discrete "time" variable j (j — 1, 2, .3, L+ 1, upper line in the table) to each 
permutation encountered in the second line and, hence, allows to order this permutation. 

Now, in a standard notation, we call j the "rise", if TTj < i^j+i, otherwise, if -Kj > ftj+i, we refer to it as the 
"descent". Further on, if we have simultaneously itj-i < ttj and itj > Ttj-i, we call j - the "peak". 

Then, the Permutation Generated Random Walk is defined by the following recursive procedure: 
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(i) at time moment j = the walker is at the origin; 

(ii) at time j > the walker makes a step to the right if j is the rise, and makes a step to the left if j is the descent. 

Note that, evidently, if j is a peak, the walker makes a left "U-turn" . 

Statistical properties of the PGRW were studied in Ref.[23|, where the distribution of the end-to-end distance and 
intermediate points of the trajectory, probability measure of different trajectories, the number of U-turns, as well as 
various correlation functions have been analyzed with respect to the uniform measure on the ensemble of random 
permutations. 

Using the methods developed in Ref. 22], and exploiting the connection between the one-dimensional ballistic 
deposition process and dynamics on permutations established in the previous section, one can straightforwardly 
reconstruct the probability distribution P{M, L) of having M peaks in the uniform ballistic deposition process in the 
planar box composed of L columns in the stationary regime. According to [Tflj |. the generating function of number of 
peaks W(s,L) can be expressed in terms of the generating function of the number of rises, U(t,L), in the ensemble 
of equally weighted permutation matrices of size L (the so-called Eulerian number), Eq.Q. The explicit relation 
between W(s, L) and U(t, L) is: 

/ 2 \ L+1 

W{ Sl L)={ — \ U(t,L) (25) 



where 



4t 2 - s - 2VT - 1 (1 - v 7 ! - s) 2 



1 (1 + t) 2 ' 
Taking into account the identity established in p^ : 



(26) 



on roo /„■ \ L+l 

A(M,L) = — cos(x(2M-L- 1)) , (27) 

7T In V X I 



we arrive at the following explicit expression for the generating function W(s, L): 

/ 9 \ L+l oo 

W(s,L) = (_±-) J2t M ^A(M,L) 



oo 



2 L + l Ll J ^ sina .^+l e -Mi+l) 4 ((l + e 2ix(3+i)) t _ e 2ix(2+i)_ e 2^ 



(28) 



n(l + t) L+1 J \ x J (1 + e 4lx )t - e 2lx (1 + t 2 ) 



with t given by (|26|l . Polynomials W(s,L) admit also another useful representation directly following from the 
relations between generating functions W(s,L) and U(t,L) found in [lflj : 



W(s, L) = 2^ 2 (1 - s)^)/ 2 LL (i+1) ( (1 ^ )2 ) 



(29) 



where Li„(x) is the Poly-log function. 

Exact expression for the distribution function P(M,L) obtained from l|29ll is given by Eq.(|SJ). More details on its 
derivation are presented in the Appendix. In the asymptotic limit L>1 the double sum JHJl converges to a Gaussian 
function with a nonzero mean see Fig El 

Calculating the first three central moments of the limiting distribution function in @ for L> 1, we get: 

L 1 
fi\ B = MP(M,L) = -L, 

M=l 6 

^ D = ((M-(M)) 2 ) = ±L, (3°) 

^ d = ((m-(m)) 3 ) = o. 
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FIG. 6: Exact distribution function Q (•) and its Gaussian approximation Q (line) for L = 90. 

Note that the first two expressions coincide with the expectation and the variance computed by J.Desbois in [Isj 
using decoupling of the hierarchy of correlation functions. The third central moment, i.e. ^3°, of the distribution 
P(M,L), naturally turns to be equal to since Gaussian limit is considered. In what follows we proceed to show 
using the method of the correlation functions (see the next Section for details) that all moments of the distribution 
may be calculated exactly, for arbitrary finite L. As a matter of fact, the value of n]P is 

^ = ((^-<M)) 3 ) = -A L , (31) 

i.e. it is not zero and grows linearly with the system size L. This value of /ig D coincides with the one computed using 
the generating function in 124|) or Q28[l. 



B. Correlation functions of the PGRW within the operator formalism 

In Ref.[23 it has been shown that a non-Markovian PGRW can be effectively described using an auxiliary recursively 
constructed Markovian process, which has the same distribution as the PGRW. Employing the ideas of Hammersley 
|26j (see also (2^ for more details), who has used such an approach in his celebrated analysis of the problem of the 
longest increasing subsequence in random permutations, the following random walk model has been constructed: 

• Suppose one has an infinite in both directions line of integers and a random walker initially at the origin; 

• Trajectory Yi of this auxiliary process is built step by step: at time moment I define a real- valued random variable 
xi+i, having a uniform distribution in [0, 1]. If xi+\ > xi the walker is moved to the right, i.e. Yi — Yi-\ + 1: 
otherwise, it goes to the left, i.e. Yi = Yj_i — 1. At the next time moment, chose compare it with xi+\ 
and move the walker accordingly, and etc. 

Two strong results have been proven in Ref. p2^ | concerning the relation between this recursively constructed Marko- 
vian process and the PGRW: First of all, it has been shown that the probability P(Yl = X) for the trajectory Yl of 
the auxiliary process to appear at site X at time moment L is equal to the probability P{Xj^ = X) for the end-point 



Exact distribution P(M,L=90). 
Gaussian approximation 
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X^ of the PGRW trajectory X^ generated by a given permutation of [L + 1] to appear at site X. This signifies 
that 

P(Y L =X)= { - ^[ L A ^ _ jL y (32) 

(Recall that A(...) is the Eulerian number). 

Second, it has been proven that the probability P{X^ = X) for the PGRW trajectory X^ at intermediate time 
moment I = 1, 2, 3, . . . , L to appear at site X is equal to the probability P(Yj = X) for the trajectory of the auxiliary 
process Yi to appear at time moment I at site X. That is, 



P{X™ =X) = P( Yl =X) = P(XP =X)= {1 - { ^ )X+l) A ( *±^, I ) , (33) 

which signifies that distribution of any intermediate point X; of the PGRW trajectory generated by permutations 
of a sequence of length L + 1 depends on I but is independent of L. 

Using this equivalence, it is thus possible to write down explicitly the probability measure of any given PGRW 
trajectory (or of some part of it) as a chain of iterated integrals over real- valued random variables xi. Recollecting the 
definition of the PGRW and noticing that each trajectory mirrors one-by-one a unique " rise-and-descent" sequence 
in a random permutation, we are thus able to represent a probability of any "rise-and-descent" sequence as a chain 
of iterated integrals over real- valued random variables x\ uniformly distributed in [0, 1]. 

More specifically, following Ref.j22|, consider a given "rise-and-descent" sequence a(L) of length L of the form: 

a(L)=( 1 2 3 • L 
\ «i «2 «3 ■■■ a L 

where ai can take either of two symbolic values: f or J.. Consequently, the first line in the table is the running index I 
which indicates position along the permutation, while the second line shows what we have at this position - a rise or 
a descent. Assign next to each symbol at position I an integral operator; for a rise (f) it is /;(T), while for a descent 
(I) it will be where 

I*(t)= I dxi and Ii{i)= I' 1 dxi. (34) 

Jxi-! JO 

To each L-step trajectory one associates next a characteristic polynomial Q(x, a(L)) defined as an "ordered" product 

M 

L 

Q(x,a(L)) = :JJ: /,(o,) • 1, (35) 

where a; = {T,l} for I = l,...,L. The statistical weight, i.e. the probability distribution function, P(a(L)), of this 
given "rise-and-descent" sequence a(L) in the ensemble of all equally likely permutations is then simply given by 

P(a(L))= [ Q(x,a(L))dx. (36) 



In may be expedient to illustrate this formal representation on a particular example. Consider, e.g., a given "rise- 
and-descent" sequence of the form {j*, j~, J., |, f}- For this sequence, the characteristic polynomial Q(x, a(5)) reads: 

1 1 x 2 1 1 

Q(x ) a(5))=/ 1 (t)/ 2 (T)4a)/4(T)/5(T)-l = J d ^J d ^J d ^J d ^J dx ^ 1 = ^-^ + ^-^ + ^ 

X Xl 0X3 x 4 

and, hence, the probability of this particular configuration is 

f 1 19 

P(a(5))= ^ Q(x,a(5))dx=— . 
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We note that construction of the probability measure in (|34|) - (|36[) can be most easily understood by considering 
the following example. Suppose there are three "markers" representing the particles with the coordinates X\, X2, £3 
(0 < {xi, X2, X3} < 1) — see Figd Markers in FigEJcan be independently deposited in the interval [0, 1] with uniform 
distribution. It is obvious that the probability P("f J.) for three particles to create a peak, is defined by the probability 
of a configuration with < x\ < X2 and x% > X3 < 1. Thus, 



P(U) = dx 3 dx 2 / d!Tl-l = -, 

JO Jx 3 JO 



what coincides with the operator expression 



(compare to l{3"l )l -!(3"5 ]l ). 



P(U)= / W)h(l)dx. 
Jo 







FIG. 7: Three markers creating a peak configuration. 



C. Exact calculation of the first three central moments of the distribution P(M,L) in Id. 



Now we apply the operator formalism for the computation of three first central moments of the probability distri- 
bution function, P(M, L), of having M peaks on a one-dimensional periodic lattice of size L: 



= (M) , 

(M 3 ) - 3 (M) (M 2 ) + 2 (M) 3 . 



ID 



(37) 



Let us introduce 

= 8{xi — Xi-\) and A s + = 
for i = 1, L, where 9{x) is the usual Heaviside step-function 

1 for x > 0, 



{Xi 



0{x) 



for x < 0. 



(38) 



Because of the periodic boundary conditions we have A x = 9{x\ — xl) and A \ 
formalism we can represent the expectation the following form 



1 1 



dxidx2-..dxL- 



(xl — xi). Using the operator 



(39) 
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In the sum above each term A i A^~ = 9(xi — x%-\)Q{xi — Xi+i) corresponds to the peak at the position Xi. All L 
terms in the sum in (|39[) are identical and independent, so it is possible to rewrite <|39[) as 



111 

L J J J 6(xi - Xi-x)6{xi - x i+1 )dx i - 1 dx i dx i+1 = L 





dxi-xdx l+1 



.o o 



dxi = —L. 
3 



(40) 



The second central moment fj^ can De calculated as follows 



/4 D 



(m 2 ) - (My 







L L 

^^A-A+A-A+^-x,)} 

i=i j=i 



1 2 

dx\dx2---dxL — -zt> 
9 



i i 



J ... J A l A+Aj A+ {5{x a - x b )} dxi-idxidx i+ idxj-idxjdxj + i - -L 



L 



r=0 



(41) 



where the notation {5(x a — Xb)} means the following: if some point x a € {xi-i, Xi, Xi+i} coincides with another point 
Xb G {xj-x, xj , xj + i}, then {S(x a — Xb)} = 1 otherwise {S(x a — Xb)} = and we sum over all configurations of diagrams 
shown in Fig[H^. For example, if in some configuration the points Xi+i and Xj coincide, then {5(x a — Xb)} = S(xi + i — Xj) 
etc. The value J^ D is the integral of the diagram r, a* is the "weight" of the corresponding diagram (i.e. the number 
of identical diagrams). The graphic representation of integrals in ea. H41|) is given in Fig|5^. We consider the system 
of length L with periodic boundary conditions. Instead of summing over L possible values of i we fix some arbitrary 
value Xi and perform averaging over all possible positions of Xj . That gives us L in front of the sum in (|41|) . The 
integral in (|41ll depends on j — i only. We enumerate all possible values of the integral by the index r, and compute 
the weight a* D of each integral J* D in the sum for j = 1,2, L. 



(a) 2-point correlations in ID 



J" 



J. 



J 



X, , X, X, 



Xj.; Xj Xj +i 



rlD x i-l %i x i+ 
J, •• • 



i« ># 

# *i m 



Gl 



(b) 3 -point correlations in ID 



m M 

Xj_ } X- 



ID %u X- X i+ 



X X s tl *j*k ^ 
ft * 0. « «E 



FIG. 8: Basic diagrams and corresponding integrals: a) Jl D for the computation of b) G\. D for the computation of /4 



The total number of diagrams is normalized: ^ a\ = L. The computed values of and of a* are collected in 



r=0 



the Table All such configurations are generated by shifting the point Xj with respect to Xi as shown in Fig[8Ji,b. 
If the shortest distance between Xi and Xj is larger than 2 (respecting the periodic boundary conditions), then the 



integrals over A i Af and A ■ A+ decouple and give J-J 13 = [-P(tl)] 2 



3" X 3 



i = i . The contribution of these integrals 



15 



r 





1 


2 


3 







l 

9 


2 

IS 


i 


ID 

a r 


2 


L-5 


2 


1 



TABLE I: Values of diagrams J r and of the corresponding weights a r for computing ji\ 



cancels by subtracting i. The number of such integrals is a\ B = L — 5. The total number of all other configurations 
is finite and does not depend on L, so the second central moment is proportional to L but no to L 2 . Substituting the 
values from the Table [I] into Q41|). we arrive at the following expression for the second central moment ji^p 



/4 D 



y ] a r 

,r=Q 



-L 



45 



(42) 



To calculate the third central moment /Zg D , 

/4 D = {{M - {M)f\ = (M 3 ) - (M) 3 - 3 (M) (<M 2 ) - (M) 2 ) , 
we proceed exactly along the same lines as above. The averaged third power of M is 



(M 3 ) 



E E E Ar A+A7 a ; a ^ A fe 
i=i j=i fe=i 



dx\dx2---d,XL- 



(43) 



(44) 



This quantity depends only on integrals over three groups of points Xj_i, a^, Xj+i, xj—i, xj, Xj+i and Xfc-i, Xfe, x/j+i 
and their mutual arrangement. To proceed, we fix some position of Xi (as it has been done for n^p) and consider the 
positions of Xj and Xk with respect to it. Using the diagrammatic approach we compute integrals G r in (|44|l for each 
three-point configuration. The corresponding diagrams G^ D for r = Q, ...,4 are shown in FiglSJa. We encounter the 
following possibilities: 

(i) The integral 144fl contains terms like Gq D = J J 9{xi — Xj)9(xj — Xi)dxidxj = 0; 

(ii) All three points Xi , Xj , Xk are separate. In such a situation the integrations over three groups of points are 
independent, giving for each group | according to (|4Tj)l . So, we have Gj D = [-P(tl)] 3 = jf- All terms of such a type 
in l|43l) cancel; 

(iii) Two groups have the common points and the third group is separated from them. In this case the separated 
integration over the third group of points gives (M). The integration over the rest pair of groups gives just the same 

result as the contribution to the second central moment. The factor 3 in front of ^(M 2 ) — (M) 2 ^ corresponds to 

three different ways (ij + k, ik + j, jk + i) to ascribe the indices to these points. Thus, the contributions from the 
three-points configurations of such a type and contribution of two-point configurations for (M 2 ) in (|43JI cancel; 

(iv) At least two pairs of groups have common points. Such configurations give a non-zero contribution. Integrals 
for such groups G]. A and corresponding weights c^P are summarized in the Table [n] where G r — J r for r = 0, 2, 3. 



r 





1 


2 


3 


4 







i 


2 


i 

■i 


17 




24 


L-37 


6 


1 


6 



r 





1 


2 


3 







i 

9 


2 

IS 


i 
■■1 




36 


18 


42 


15 



TABLE II: Values of integrals G* D and of the weights c* D (left) for the three-point configurations and values of integrals J* 
and weights fo* D (right) for computing [i™ . 



We can split nontrivial (r ^ 1) three-point configuration (c£ ) into the two-point configuration by deleting of 

4 3 

the one-point group in three different ways. Thus, 3 ^2 c r D = S = l^i where &J D is the number of the 

9=0 



r=0 



two-point configurations of type q obtained as a result of splitting of all possible three-points configurations. The 
two-point configurations are enumerated in the Table H*H The contribution of the two-point configurations are 
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not compensated by an appropriate term from the three-points configurations, so we have to take it into account 
manually. So, substituting expansion of (|44|l into l|43[) . we obtain 



M p = (M 3 ) - (Af) 3 - 3 (AO ((A/ 2 ) - (A/) 2 ) 



EM 

r=0 
2 . 



1D -L 
27 



q=0 



945 



-0.0021164L. 



(45) 



It is worth mentioning that the moments in Eas. (|40(l . l|42|l and 1|45JI . calculated using the operator technique, 
coincide with the ones obtained on the basis of exact combinatorial approach (|23|l . 



V. CONDITIONAL PROBABILITY p(l) OF TWO PEAKS SEPARATED BY DISTANCE I 



We aim now at evaluating the probability p(l) of having two peaks separated by a distance I, under the condition 
that there are no peaks (i.e. sequences | 1) on the interval between these peaks. According to [22J, this probability is 
given by 

p(l)= [ dxY^Q{x), (46) 
Jo 

where the sum is taken over all possible peak-avoiding rise-and-descent patterns of length I between two peaks, while 
Q(x) denote the Q-polynomials corresponding to each given configuration (see the explanations above). 

There are several possible peak-avoiding "rise-and-descent" sequences contributing to such a probability. These 
sequences are depicted in FigEl The first peak is located at position, the second peak is located at I position. We 
have fixed a descent at 1st position (variable Y) to keep a peak at the Os position (variable X). 




J I I I L 



(b) H 

X Y x } x 2 



J I I L 



J L 



( c ) V- 1 M I I J I I I I XL LX 

FIG. 9: Rise-and-descent patterns contributing to the conditional probability of having two closest peaks at distance I apart 
from each other. Configuration (a) has only one descent between two peaks. Configuration (b) has two descents following the 
first peak, i.e. a "through" at xi, and (c) presents a generalization of (b) over configurations having m descents, (m = 1,2, ... ,1), 
after the first peak which are followed by I — m rises, i.e. a "through" at x m -i, 
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Now, the Q-polynomial associated with the sequence (a) in Fig|5]has the following form: 



pi pX pi pi pi pXi-i 

Q a (x) = dX dY dx x I dx 2 ... dxi-i / dx u 

Jx JO JY Jx\ Jxi-2 JO 



(47) 



Performing the integration over in the latter expression, we denote the iterated integrals over the variables Xk, 
k = 1, 2, 3 . . . , I as 

Mi{Y)= ( dx x \ dx 2 ... f xidxi. (48) 
Jy Jxi Jxi-i 

Notice now that Mi obey the following recursion scheme: 

M ; (Y ) = J M ; _i (X)dX, M {Y)=Y. (49) 
Introducing the generating function of the form: 

oo 

M(r)=^M,(yy, (50) 

we readily find that it obeys 

M(Y) - Y = z { M(X)dX, (51) 



and consequently, Mi are simply the coefficients in the expansion 



oo 

M{Y) = Miz 1 = - 1 - (1 - z) exp(z(l - x)) , (52) 



1=0 

which are given explicitly by 



(~\ — Y\ l C1 - YY+ 1 1 



Hence, the Q-polynomial associated with the configuration (a) in Fig|5] obeys 

Q a {x) = [ dX [ M l _ 1 (Y)dY = i f dX [ [(1 - Y) l -\l -l + Y)] dY, 

Jx Jo Jx Jo 

and the contribution of this very configuration to the probability p(l) reads: 



(54) 



P a {l) = y j dx ( dX [ [(i-y)'-i(z_i + y)] dY, (55) 

Jo Jx Jo 

Next, we turn to the contribution coming out of the general configuration (c) in FigEl The Q-polynomial associated 
with this configuration of rises and descents for a fixed m is given by 

,1 rX rY r-Xl rX 2 j-X m -2 fl fl r%l-\ 

Q c (x,m) = / dX / dY / dx\ \ dx% I dx$ . . . / dx m -i / dx m . . . / dxi-i / dx; 

Ja; JO Jo Jo Jo Jo Jx m -i Jxi-2 Jo 

pX pY pX\ pX 2 pX m —2 pi pi 

dX / dY / dx\ I dx 2 I dx% . . . dx m -i / dx m . . . xi-idxi-i 

Jo Jo Jo Jo Jo Jx m -i Jxi-2 

fX j-Y r-Xl rX 2 rX m -2 

dX I dY I dxi / dx 2 / dx 3 . . . / Mi- m (x m -i)dx m -i. (56) 

Jo Jo Jo Jo Jo 

Let us note, that Q a (x) = Q c (x,m = 1). Consider next a recursion scheme of the form 

N m (Y) = z [ N m _ x {X)dX, (57) 



where Nq(Y) is some arbitrary function <&(Y). Introducing the generating function 

oo 

N(Y) = N m (Y)z m , 

we get that it obeys 

N(Y) - $(y) = z ( N{X)dX. 



Solution of the latter equation can be readily obtained by standard means and reads 



N(Y) = $(0)exp(zy) + 



dX 



exp 



(z{Y-X)) 



dX. 



Finally, expanding the rhs of the latter equation in powers of z, we get that N m (Y) are given explicitly by 



N m (Y) = $(0)— + 



Y m f Y d®(X) (Y - X) 



dX 



dX = / $(X) 



(Y-X 



i rn — 1 



(m-l)! 



-dX. 



Now, we notice that, as a matter of fact, the multiple integral over the variables Y and x^, k = 1, 2, . . . ,m 
the rhs of Ea. l|5f>J) becomes just the function N m (Y), if one takes $>(Y) = Mi- m (Y). This implies that 



A 



dY / dx\ I dx2- 
a Jo 



Xi 



%m-2 



Mi- m (x m -i)dx 

r, 



X 



(X - F) m " 1 (l - Yy- m (l -m + Y) 
[m- 1)!(Z- OT + l)! 



dY. 



x 



Substituting to — 1 we obtain here J Mi_\dY . Consequently, we find that the desired probability obeys 



p(l) = dx dX 

JO J x 



X 



£ 

m — 1 



{X - F) m_1 (l - Yy- m (l -m + Y) 



(m - — 771 + 1)! 

Let us introduce the generating function F(z) — z 'p(0- We can represent the sum 



dY. 



1=2 



oo oo 



]T E zl f( l > ™) = E E z h+m f(k + m, to), 



1=2 m=l 



m= 1 k= 1 



where k = I — m. Using 



and 



E 



(x-y) m - 
(m- 1)! 



_ 7p z{X-Y) 



|«'<» + ^^r-*- Y '( 1 -;) + i- ir - 



we obtain 



oo / -i 

*•(*) = £**(/ 

Z=2 V- 70 



/ dX 
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Hence, the generating function is given by 

Note that the numerical values ^ (for I — 2) and i (for I — 3) correspond to the values obtained in j23|. In general, 
we get the following explicit expression 

P (l) i - 2-^- + J^J) ^ l -^ l + 2) (69) 

P[L) ~ 2\(l + l)\ \l + 2)\ + (l + 3)\J 1 (1 + 3)1 ' (bJj 

which defines the probability p(l) of finding two peaks separated by the distance I with no peaks between these two 
points for arbitrary I. The function p(l) is depicted in Fig llUI 



. 14r 

0.12 
0.1- 
0.08- 
0.06- 
0.04- 
. 02 - 




FIG. 10: Probability of having two peaks being at distance / apart of each other with no peaks inbetween. It has a maximum 
at I = 2. 

Note also that p(l) in ea. (|69|) coincides with the distribution function of the distance between two "weak" bonds 
obtained by Derrida and Gardner pl| in their analysis of the number of metastable states in a one-dimensional 
random Ising chain at zero's temperature. 



VI. PROBABILITY DISTRIBUTION FUNCTION P{M,L) IN 2D. 



In this last section we are going to study the statistics of peaks in a two-dimensional BD landscape h(i,j,n) by 
converting this problem to the analysis of the distribution of peaks in the associated permutation matrix. 

It is convenient to represent the two-dimensional base (i, j), (i, j = 1, 2, 3, . . . , m, m x m = L) as a one-dimensional 
set with long-ranged correlations. Namely, reading the lattice from the left to the right in the line and line-by- 
line from top-to-bottom, exactly as an electron beam does to highlight the image on the TV screen, we can rewrite 
equations Q and J3J| renaming h(i,j,n) as: h(i, j,n) = h(k,ri), where k = 1,2,3, . . . , L. Hence, h(i — l,j',n) = 
h(k — 1, n)\ h(i + n) = h(k + 1, n); h(i, j — 1, n) = h(k — m, n); h(i,j + 1, n) = h(k + m, n), where 1 < k < L. 

Given Eq. @ , we can describe the growth of a two-dimensional BD landscape over the base m x m as a stationary 
"updating dynamics" in the L x L permutation matrix with the uniform distribution on updating events. The 
uniform "updating dynamics" on permutations generates the PDF of peaks identical to the PDF of peaks computed 
over the ensemble of all L\ equally weighted permutations (compare to the one-dimensional case). Thus, repeating the 
arguments of the previous section, we can construct an analogue a two-dimensional PGRW for permutation matrix 
with finite-length correlations. 

The most significant and the only difference between one- and two-dimensional models is in the definition of a 
" peak" . In ID case, a peak appears at the position k of the permutation matrix if the corresponding permutation iik 
is larger then the nearest neighboring permutations itu-i and 7Tfe-|_i. In 2D case, k is a peak if and only if tt^ is larger 
than 7Tfe_i, 7Tfc+i, TTk-m and Tik+m simultaneously. In the next subsection we generalize the operator formalism over 
the 2D case and calculate the PDF of peaks for the uniform BD process in two dimensions. 
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A. Central moments of the probability distribution function in 2D. 



To obtain the limiting (L — ► oo) Probability Distribution Function P(M, L) of having exactly M peaks on a two- 
dimensional square base L = m X m, we calculate the first three central moments of the distribution function and 
then construct the corresponding Edgeworth series (cumulant expansion) |23|. This enables us: a) to show that in the 
limit L — ► oo the function P(M, L) converges to the Gaussian distribution, and b) to present an explicit expression 
for P(M, L) in this limit. 

The total number of points for 2D system is L = m 2 . Applying to the 2D case the same arguments as in ID case, 
we define 



0(xi,j — Xi—ij) for the "northern" neighbor 



*0 = e ( x i,i 
= 6(xi t j 



A 1 ■ 



Xi j+i) for the "western" neighbor 
Xi+tj) for the "southern" neighbor 
Xij-i) for the "eastern" neighbor, 



(70) 



(compare J7D|| to J3S)). 

Using the operator formalism we obtain for the first central moment 
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(71) 
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We define dx 



As in ID, all terms in l|71|l are identical and independent. 

2. The same method is used for computing the second central moment. 
dxi^ijdxij+idxi+ijdxij-idxij for integration over the point x^j and its 4 neighbors to make the expressions more 
compact. Instead of summing over we fix the position of the first point Xij and perform the summation over 
all possible positions of x k ^ with respect to Xij. Then we enumerate all different integrals J 2B and compute the 
corresponding weights a 2B . The configurations which contribute to J 2B are shown in Figlllb.b. As in ID case, the 
(5-functions 5{x a ^ — x c ,d) cu t off the diagrams with coinciding points (see the explanation after 141|l ). 

The values of integrals J 2U and their weights a 2B are collected in Table ITO! 
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TABLE III: Values of integrals J^ D and of the weights a^P for computing [pp . 



The total number of integrals is a r° = L = m 2 . The contribution from the integral Jf D is exactly canceled by 
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the term — i coming from two independent points: 
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(a) 2-point correlations in 2D 



(b) 3-point correlations in 2D 
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FIG. 11: Basic diagrams and corresponding integrals: a) J^ D for the computation of ^2°; b) for the computation of /i| D . 

3. For computation of the third moment in 2D we use again the same method as in ID. Namely, we fix the first 
point Xi y j and enumerate all possible configurations of three points. Different types of three-points configurations 
are depicted in Figlllb). Other configurations with the same contribution and the same topologies but slightly 
different conformations are not shown in the figure. The integrals G^ D and J% have the same values for r = 0, 2, 3, 4. 
In the configuration Gf^ all three points coincide, so Gf 3 = jL The most general form of the integral G^ D is as 
follows 

l l 

= j ■■■ J A^.A^.A^.A^ A^A^A^A^ Aj^A^A^A^ {S(x a , b - x c , d )} rf Xiji dx fe ,^x m ,„, (73) 



where e?Xij = dxi-\jdxi,j-\dxi+i,jdxi t j+\dxi.j and the (S-functions in H73J1 cut off the coinciding points. For example, 
if in some configuration the points Xi_j + \ and Xk.i-i coincide, then we include the function 5(xk,i-i — Xi_j + i) etc. 

It is possible to simplify integrals of such a type by changing the limits of integration. For example, 
11 l 

f f Aj jdxijdxi-i.j = J Xijdxij . In such a way the integral G§ D can be expressed as follows: 
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(74) 



The values of integrals Gq D — Gg D are collected in the Table llVl 

The symmetry of configurations defines the weights c r . For example, the diagram of the integral G| D in Figlllb 
has two orientations along vertical and horizontal lines. Thus the total weig ht of the integral G| D is cf = 36. 
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TABLE IV: Three-point configurations for computing ^3°: integrals G> D and weights cJP (left) and integrals J^ D and weights 
b™ (right). 



The values of c r are given in the Table HVI Totally there are c 2 , = 313 different nontrivial configurations. We 
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r/1 



can split each of the nontrivial three-points configuration as it has been done in ID case. There are three ways to 

do it. Values of b r of the corresponding two-point configurations are shown in the Table Hvl with the total number 
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VJ 6 2D _ ggg £ g^fa configurations. Now it is possible to compute the third moment 
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B. The distribution function P(M,L) in the limit L — > 00. 



The second and the third cumulants of P(M, L) in 2D are equal to the second and the third central moments: 

(76) 



Introducing the normalized deviation, x = — 
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P [Hi + xa, Lj in a form of the Edgeworth series (cumulant expansion) [2S 



p(x, L) ~ 3 (x) ( 1 + -L/C*) + o ( _L ] ) , ( 77 ) 



where g(x) is the Gaussian function g(x) = eT x / 2 , and f(x) is defined by /i 2 , and /i| D (see for details) 
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Substituting in (JTSJ), we get 



fix) = VZ— «-(*»- to) (78) 

( K 2 ) 



/W — f^V 72 ^- 3 -^). (79) 



32175 V 13 I 



Let us emphasize that f(x) does not depend on L, as one may readily check by substituting exact results obtained 
for the cumulants to Ea. (|78|l . To verify the expression in (|79|l . we have performed the numerical simulations for the 
discrete 2D permutation-generated model with the periodic boundary conditions and have computed the distribution 
function p(x, L) numerically. In Fig 1 12b we present the data of the numerical simulations for p{x, L) and compare it 
against the Gaussian function g(x) for system sizes L — 100 and 256. Furthermore, in Figll2b we plot the ratio ^ferp 
as the function of x, which shows that the deviation of the numerically computed function p(x, L) from the Gaussian 
function g{x) is actually very small. 

One clearly sees that for L — > 00 the normalized probability distribution function p(x, L) converges to the Gaussian 
function g{x). 
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FIG. 12: (a) Results of numerical simulation of p(x,L) for L = 100, 256: (a) comparison of p(x,L) with the Gaussian function 
g(x); (b) comparison of the function p(x,L)/g(x) with 1 + -jff(x). 



VII. CONCLUSION 

To conclude, we have studied the probability distribution function P(M, L) of the number of local peaks in one- 
and two-dimensional surfaces obtained in a standard model of next-nearest-neighbor ballistic deposition process. Our 
analysis was based on two central results: 

(i) A proof, presented in this paper, of the fact that a ballistic deposition process in the steady state can be 
formulated exactly in terms of "rise-and-descent" patterns in the ensemble of random permutation matrices, which 
made it possible to interpret the BD model in terms of permutations of the set 1,2,3, L, where L is the number 
of lattice sites and the numbers drawn at random from this set determine local heights of the surface. 

(ii) An observation made in Ref. [2^| that the "rise-and-descent" patterns can be treated efficiently using a recently 
proposed algorithm of a Permutation Generated Random Walk. 

In one-dimensional case we have found a closed-form expression (|28|l for the generating series Wis, L) — 
LI Em=o s m+1 P(M, L) of the distribution function P{M, L) of peaks in a bounding box of size L. Inverting this 
expression, we got an exact and asymptotic forms of P(M,L). 

Besides, in one-dimension we calculated the probability p(l) of having two peaks separated by a distance I, under 
the condition that there are no peaks in the interval between these two peaks. The function p(l) is given by expression 

jug. 

For two-dimensional case, we reformulated the BD process in terms of an "updating dynamics" on permutations 
with certain finite-range correlations. Using this approach, we extended the operator formalism of p2l | to such 
correlated permutations and evaluated three first central moments of the PDF. Then, we have obtained P(M, L) in 
the asymptotic limit L — > oo using expansions in the Edgeworth series [23| — see (|77|I - H79|) . We have shown that in 
2D the PDF also converges to a Gaussian function as L — > oo. 
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APPENDIX A 



Let us exploit the relation between the generating function of peaks, 

oo 

W(s,L)=L\ s M+1 P{M,L) 



M=0 



and Poly-log function (see [lij for details): 
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Introducing s = (i^a > rewrite 1|A1|) taking into account that t = 
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Use now the expansions 
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and rewrite W(s, L) in 1A2JI as follows 
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After resummation of (|A5|) we arrive at the following expression: 
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Hence, by definition of the generating series, we get 
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what gives the desired final expression for PDF P(M, L) in one dimension. 
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